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We propose a model of inflation driven by minimal extension of SUSY, commonly known as MSSM. 
Starting from gauge invariant flat directions in the n = 4 level comprising of QQQL,QuQd,QuLe 
and uude, we construct the inflaton potential and employ it to investigate for its consequences 
around the saddle point arising from the non-vanishing fourth derivative of the original potential. 
To this end, we derive the expressions for the important parameters in MSSM inflation using the loop 
corrected potential. We further estimate the observable parameters and flnd them to flt well with 
recent observational data from WMAP7 by using the code CAMB. We also explore the possibility 
of primordial black hole formation from our model. Finally, we analyze one loop RGE and compute 
different phenomenological parameters which could be precisely determined in LHC or future Linear 
Colliders. 



I. INTRODUCTION 

The paradigm of primordial inflation is, by far, the 
most satisfactory explanation for early universe phenom- 
ena i- As a general prescription, inflation occurs due 
to a slowly rolling scalar field, the inflaton, dynami- 
cally giving rise to an epoch of accelerated expansion 
dominated by a false vacuum 0. Primordial quantum 
fluctuations of inflaton are responsible for creation of 
matter content and observed perturbations in the Cos- 
mic Microwave Background Radiation (CMBR). Further, 
slow-roll inflationary scenario generically predicts almost 
Gaussian adiabatic perturbations with a nearly flat spec- 
trum, which conforms well with the latest observations. 

Recently, some interesting proposition of inflationary 
model building was brought forth by Minimally Super- 
symmetric Standard Model (MSSM) where the inflaton is 
a gauge invariant [3, 4J n = 4 level combination of scalar 
superpartners squark and slepton fields and fermionic 
superpartner gauginos which are candidate Cold Dark 
Matter (CDM) particles. However the original poten- 
tial for n = 4 level is unable to extract a suitable sym- 
metry along the fiat direction. To serve this purpose 
the usual way is to incorporate saddle point mechanism 
to the MSSM potential leading to vanishing of the sec- 
ond derivative and the slow roll phase is driven by the 
next leading order derivative of the potential [l|, [^-[ll ■ In 
most of the phenomenological situations, a fine tuning 
mechanism is needed to place the flat direction field to 
the immediate neighborhood of the saddle point. It is 
worthwhile to mention that MSSM inflation occurs at a 
comparative lower scale. This is in strong contrast with 
the conventional class of models where the unfamiliar in- 
flaton couplings to Standard Model (SM) are originated 
through arbitrary gauge singlets leading to the field mag- 
nitudes at GUT scale or higher, and hence, face problems 
in satisfactory quantitative estimation of a huge sector of 
the post-inflationary evolution i.e. thermal history of the 
early Universe, baryon asymmetry and CDM. Herein lies 
the most appealing feature of MSSM inflation for which 
known SM couplings are measurable in laboratory exper- 
iments such as Large Hadron Collider (LHC) ff\ or future 



linear colliders. 

In the present article we will consider a specific MSSM 
scenario where, for a specific choice of soft supersymme- 
try (SUSY) breaking parameters A (trilinear couplings) 
and the infiaton mass m^, the potential is D-fiat along 
the QQQL,QuQd,QuLe and uude directions. For our 
model existence of saddle point is guaranteed by the non- 
vanishing fourth derivative of the potential, which makes 
the potential more fiat than the previous ones. This im- 
plies more precise information in the RG flow. As we 
will show, this is the highest level of precision constraint 
one can impose on RG fiow keeping the effective poten- 
tial renormalizable in the vicinity of the saddle point. 
Our primary intention is to investigate for the analytical 
as well as the numerical expressions for different obser- 
vational parameters for MSSM inflation with these new 
flat directions. As it will turn out, they match quite well 
with latest observational data from WMAP7 Q and are 
expected to fit well with upcoming data from PLANCK 
Q . Additionally we have explicitly shown the connection 
between running and running of the running of spectral 
index to the Primordial Black Hole (PBH) formation. To 
this end we get the fine tuned parameter space which is 
also in good agreement with present estimates of cosmo- 
logical frameworks. We have further explored features 
of the MSSM from the solution of one loop RGE which 
could be measured by LHC or future linear collider. 



II. FLAT DIRECTIONS AND POTENTIAL 
AROUND SADDLE POINT 

Let us start with n = 4 level superpotential [l^ 
= fE?ti a/(QQQL)/ + /3/(QuQd),- 



E/=i 7/(QuLe)/ + Y.i=i (5/ (uude)/ 

(1) 

The renormalizable flat directions of the MSSM at n=4 
level correspond to the gauge invariant monomials sub- 
ject to the four additional complex constraints [l^ two 
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each from 



= ^^H2 + Xtj^Q^Ub = 0, (2) 



-^iH: + V^QZd, + X'iL^eb = 0, (3) 



which can Hft the flat directions which do not contain a 
Higgs field. Here Xu,Xd and are the Yukawa cou- 
plings, Hu,Hd are the Higgs superfield and the /x- term 
appears in the renormalizable part of the superpotential 
of MSSM . Consequently the equation([T|) breaks into four 
parts, each one of them now being flat: 
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resulting in W^]'' « y^j-^ (1,2,3,4)). Consider- 

ing any one of the above flat directions leads to the one 
loop corrected effective potential 
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and in Gmssm = SU(3)c «> SU(2)l «) U(1)y the rep- 
resentative flat direction field content is given by 
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(9) 

Here mo,Ao and A4_o are the values of the respec- 
tive parameters at the scale /io and Di , D2 and D3 
{\Di\ ^ IVi) are the fine tuning parameters. Addi- 
tionally in the field contents l<Bi,B2,B3<3 are 
color indices, 1 < a, b, c < 3 denote the indices for 
quark and lepton families and 1 < Ii, I2, 13, 14 < 2 are 
the weak isospin indices. The flatness constraints re- 
quire that Bi 7^ B2 7^ B3 for quarks, Ii 7^ I2 7^ I3 7^ I4, 
Ed=i Pd 1 e K for leptons and a 7^ b 7^ c for 
both. In eqn([8]) represents the soft SUSY breaking 
mass term, the radial coordinate of the complex scalar 
field $ = (j)exp{i9) (g C) and the second term is the so 
called A-term which has a periodicity of 27r in 2 D along 



with an extra phase 9a- The radiative correction slightly 
affects the soft term and the value of the saddle point. 

For n = 4 we get an extremum for the principal values 
of 6* at 61 = (JUI^ (where m € Z) 
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(10) 



which appears from the constraint V {(po) = as a nec- 
essary condition for saddle point. However, this condition 
alone will not lead to saddle point. Rather, we have to 
make the potential sufficiently flat which can be achieved 
by vanishing higher derivatives of the potential. In this 
article, we consider non-vanishing fourth derivative of the 
potential resulting in saddle point. This will imply more 
fine-tuning but increased precision level in the informa- 
tion obtained from RG flow. Below we demonstrate how 
this is materialized. 

As discussed, V (0o) < will give us secondary local 
minimum. This leads to constraint relations: 
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and V"'{(j)o) = 0. In 
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_(3+-D3) 

G2 = [{1 + Di) {3 + - (1 + 3Di) {1 

-1 



G 
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limit l^il < l,|£i2| < 1 and {Ds 
(fit"" [^ + ^~^Y and A - 
where 
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For the 



itree 
"0 



^ 1 which gives 0o = 

^tree L"^ ^ 2 6 J ' 

^^i^^ and Atree = 2^m^{ct)o) rep- 
resents tree level expressions. This means, during RG 
flow mentioning two parameters only {Di and D2) will 
suffice instead of the usual three parameters in earlier 
MSSM models. This results in more precise informa- 
tion in RG flow. One may get tempted to vanish further 
higher derivatives of the potential in order to evaluate 
other unknown parameters (Di and D2) without going 
into RG flow but this will make the effective inflaton po- 
tential in the vicinity of saddle point non-renormalizable. 
So, this is the highest level of precision constraint one 
can impose on RG flow parameters. 

Consequently, around the saddle point (f>Q , the inflaton 
potential can be expanded in a Taylor series as. 



V{^) = Co + (74 



(13) 
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where Cq = V{(j)o 
l + Alog(|)" 
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In what follows we shall model MSSM inflation with 
the above potential. 



III. MODELING MSSM INFLATION AND 
PARAMETER ESTIMATION 
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FIG. 1. Variation of equation of state parameter(tj(a;)) versus 
shifted inflaton field (x) 



For brevity, let us introduce a change of parameter (jj 
X — (f) — (j)Q which represents the inflaton with shifted 
origin. Using this new notation of field the slow roll pa- 
rameters are given by, 
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(14) 
(15) 
(16) 
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where a prime denotes d/d4) = d/dx. During slow-roll 
inflation £„, |77u|, \£^y\, |cr«| ^ 1 and the end of the inflation 

corresponds to \xf \ ~ ( sm"c^ J ''^herex/ — (pj — tpo- In 
this context equation of state parameter can be expressed 



-C^^x^ + iiClM^x<'-2CoCiX*-C^, 



Clx» + ^ClM^x<i+2CoCix'^+Cl 



(18) 



which implies the energy scale of MSSM inflation is 

- (0.409 - 1.301) X lO-^M expUcitly shown in the 
allowed region in Fig(IT]). 

The number of e-foldings are defined [ll| for our model 
in the limit Co ^ C4 as 



TV log f 1 f^' 
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Further, in this framework the expressions for ampli- 
tude of the scalar perturbation, tensor perturbation and 
tensor to scalar ratio are given by 
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Here a;^ represents the value of the inflaton field at the 
horizon crossing. For our model expression for spectral 
index, running and running of the running reduces to the 
following form: 
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FIG. 2. Variation of the scalar power spectrum(As) vs scalar 
spectral index(ns) 
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CMB TT Angular Power Spectrum 



FIG. 3. Parametric plot of the logarithmic scaled ampli- 
tude of the scalar fluctuation (Zn(As)) vs logarithmic scaled 
amplitude of the running of the spectral index (ln(|Qs|)) . 



Fig© depicts the behavior of the scalar power spec- 
trum as a function of scalar spectral index. For M = 
70 the scalar spectral index is within the bounds of 
WMAP7+BAO+hO data for model kCDM+sz+lens 0. 

FigQ shows the behavior of amplitude of scalar fluc- 
tuation as a function of running of the spectral in- 
dex. For the best fit value of Co = 2.867 x IQ-'^^M^, 
Ci = -1.685 X 10"^^ and Af = 70 the cosmological pa- 
rameters obtained from our model is = 2.498 x 10~^, 
A? = 1.936 X 10-39, Hs = 0.957, iit = -1.550 x lO'^o, 
r = 1.240x10-29, as = -0.612x10"^^ Ks = 1.749xl0"-l 

Further, we use the publicly available code CAMB 
to verify our results directly with observation. To 
operate CAMB at the pivot scale fco = 0.002 Mpc"^ 
the values of the initial parameter space are taken for 
Co = 2.867 X lO-^^Af* and TV = 70. Additionally 
WMAP7 years dataset 8] for ACDM background has 
been used in CAMB to obtain CMB angular power spec- 
trum. In TablcUwe have given all the input parameters 
for CAMB. TablcHI] shows the CAMB output, which is 
in good agreement with WMAP seven years data. In 
FigHwe have plotted CAMB output of CMB TT angu- 
lar power spectrum Cf^ for best fit with WMAP seven 
years data for scalar mode, which explicitly show the 
agreement of our model with WMAP7 dataset. 
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TABLE I. Input parameters 
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TABLE II. Output obtained from CAMB 




FIG. 4. Variation of CMB angular power spectrum C^^ for 
best fit and WMAP seven years data with the multipoles 1 for 
scalar mode 



Now in the context of any running mass model one can 
expand the spectral index with the following parameter- 
ization Il3l: 



r^(7^) = n^(fco)- 



0:2(^0) 



ln(fco7^)- 



In^ (koTZ)- 
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(27) 

with TZ <C 1/fco, i.e. \n{koTZ} < 0. This is identified to 
be the significant contribution to the Primordial Black 
Hole (PBH) formation. Here the parameterization index 
z : [s{scalar),t{tensor)] and the explicit form of the first 
term in the above expansion is given by 



I nt{ko) li z = t. 



(28) 



Existence of the running and running of the running is 
the key feature in the formation of PBH in the radiation 
dominated era just after inflation 14] which could form 
CDM in the Universe. The initial PBHs mass A^pbh 
is related to the particle horizon mass M. by TMpbh = 
= ^jpn-^ at horizon entry, TZ = [oH)-^ . This is 
formed when the density fluctuation exceeds the thresh- 
old for PBH formation given as in Press-Schechter theory 
by Bm 



f{>M)^2j der{0;Min)). (29) 

Here V{0:M(TZ)) is the Gaussian probability dis- 
tribution function of the linearized density field 
smoothed on a comoving scale TZ by T'{0;TZ) — 

\/2ttS (u) "^"^P (~ 2i;f (R) ) '^^cre the standard deviation 
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For our model power spectrum for ry) is given by 



8JV3CoGfc'^))'^(l + fc^i)^) 



where Ks('7) = 



(31) 

is the 



equivalent expression for running of the running in terms 
of conformal time 77 and J = -j^-pr^^^- Additionally 

and G = Substitut- 



we have used $(/) 



^3Cn 



ing ean ipTj) in eqn([5ni) and using eqn (P5| at the hori- 



zon crossing we get S0{TZ) 



25Af7r2Ks ^ I- 2 

for UsiTZ) > 3. Consequently eqn (P^ gives /(> A^) 



8V3Co.7GT^r("sCR.)+3)i 



7erf 



Here we fix 7 ~ 0.2 during the radi- 



ation dominated era [16l| for proper numerical estima- 
tions. In general the mass of PBHs is expected to depend 
on the amplitude, size and shape of the perturbations 
p7| . As a consequence the PBH mass is given by [l^ 

A^PBH = iMcqikoq'R')'^ (^gf^) ^ ' ■^here the subscript 
"eq" refers to the matter-radiation equality. Here we 
use 5* = 228.75 (all degrees of freedom in MSSM), while 
.9*,oq = 3.36 and k^q = 0.0717„i/i^ Mpc~^ (Here we use 
^ 0.2670 from the CAMB output). Consequently 
the relation between comoving scale and the PBH mass 
in the context of MSSM is given by 



^ 1.250 X 10-23 ^-^PBH 



1 Mpc 



V 1 



(32) 
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FIG. 5. Variation of the fraction of the energy den- 
sity of tlie universe collapsing into PBHs as a function of 
the PBH mass, for three different values of the threshold 
0th = 0.3{dashed),0.5{solid),0.7{dotdashed). 

Fig(IS|) shows the behavior of Press-Schechter mass 
function with respect to PBH mass. With the values 
of the parameters as obtained earlier, we have M pbh — 
IQ^^gm and the corresponding fractional energy density 
/ = 0.170. 



Finally, the reheating temperature for our model turns 

( 



out to be T^h 



30a;: 



For TV ^ 
2.114 X lO^^GeV^ which is 



12007r2C|A2. 



70 it is estimated as Trh 
obviously significant input to choose the fine tuned initial 
conditions for RGE flow discussed in the next section. 



IV. ONE LOOP RG FLOW 

For the flat direction QQQL,QuQd,QuLe,uude the 

soft SUSY breaking masses can be expressed as 



(™^)qqql = i(m|^ + ml^ + ml^ + mlj 
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where 1 < a, b, c < 3 and a 7^ b 7^ c. After neglecting 
the contribution from the all Yukawa couplings except 
from the top we can express the one-loop beta func- 



tion as 
1 v^3 



1 



Z]a=l5al"^aPXQa whcrC Xq 



2^ / 'o,—i Hn\"''a\ -'^aa wiiv^xv^ -i^aa ^r6 the quadratic 
Casimir Group Invariants for the superfield defined in 
terms of Lie Algebra generators T*^ by (T"T");j = :^aaS^ 

and A = /^|l- 




FIG. 6. Running of gaugino mass {mi{fi)) in one loop RGE 
for MSSM with the logarithmic scale logjg (/i). Here we have 
used fio = 2.6 x lO'^GeV , mi(/io) = 7.546 x IQ-^TeV, ( = 1 
V i. 



In the context of MSSM 



3Y 



(for each with weak hyper charge Y^), 



X2a = | (for *„ = Q,_L,H„,Hd), 

= (for *a = u, d, el, 
Xsa = I (for *a = Q,u,d), 

= 0(for =L,e,Hu,Hd), 
where Xiq, X2a and Xsq are applicable for 
U(1)y,SU(2)l and SU(3)c respectively. So for 
the flat direction content QQQL,QuQd,QuLe,uude 
we have the following beta functions: 
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10 15 

FIG. 7. Runnine; of soft mass squared ratio ( — 4i — ^ 

one loop RGE for MSSM with the logarithmic scale logig (m) 
where /xo = 2.6 x IffGeV for C = 0.5, 1,2 for n = 4 level 
QuQd V i. Similar plots can be obtained for QuLe, QQQL 
and uude flat directions also V i. 



(a) For Soft mass 

AK)qqql = ^ {iml^ + ml^ + \Af?) 
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/*(™rf,)uude — 2^ 



^3 



) {^h'Y 



(b) For Trilinear A- term: 

■Aab _ 3(l+.5„3)A?f /'\33^2 



|5il-3p) 



IB - -4?? Vf5^l"^lr + i3ll™2p) , 
(c) For Fourth level Yukawa coupling: 

,;\aa _ 3(l+^„3) ( \33\^ >-u" ( 13 „2 , 3 ^2 
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§51) 



■ xaa _ ^ (3 2 , 3 2\ 
t^-^E — \292 + 293) 

where all the superscript a and b represent generation or 
family indices run from 1 to 3 physically representing the 
first, second and third generation respectively. For the 
one-loop renormalization of gauge couplings and gaugino 
masses, one has in general 



n • 9a^^a 



(34) 



where Xqg quadratic Casimir invariant of the group [ 
for U(l) and N for SU(N) ], I„a is the Dynkin in- 
dex of the chiral supermultiplet $a [ normalized to ^ 
for each fundamental representation of SU(N) and to 
3Y^/5 for U(1)y ]■ For the above mentioned flat direc- 
tion the running of gauge couplings and gaugino 



masses (mi(^)) obey [l^, figi = ^g,f , /i ( g^) — ^ ^ * 

where for i = 1(U(1)y), 2(SU(2)l), 3(SU(3)c) here 
di = g^j'^2 = ■g^,d3 — ^8§5- which is the simpler 
version of the equation ([34]). Now to show explicitly 
that the contributions from the top Yukawa coupling 



(A^^) are very small for an induced electroweak group 
G_EW=SU(2)^ (g) U ( 1 W breakdown, let us start with 
the Higgs potential 



Vmggsi'a,'il)=ml\U\^ + ml\II\' 



I (HH + HtHt 



(35) 



where H = Hu and H = Hd represent the Higgs super- 
fields and the relative vev of the two Higgses are given 

by 



(36) 




gi) 008(29) 



with tan{9) ~ Here 9 represents an angular param- 
eter which parameterizes MSSM. For the sake of conve- 
nience let us now write cos{29) appearing in equation([36]) 
introducing new parameterization as j20| cos(29) — ™ 



where 



<H> 



w 



/l-(<^) 

coupling can be expressed as A 
f and the top mass 43 GeV< 



Consequently the top Yukawa 
where < 9 < 



33 _ 
u — 



vsin{0) 

_ mu < 170 GeV< hgut 
comes from the RG fl ow fisl . It is evident from the above 
parameterization [20l422| that as w — ?> 1, — >■ ^ which 
implies (H) and (H) is very large and have the same or- 
der of magnitude. As a result the relative vev v is also 
large and the top Yukawa coupling is very very small for 
which one can easily neglect it from the RG flow at the 
energy scale of MSSM inflation as mentioned earlier. The 
consequence of the large vev of Higgs field can be taken 
care of by introducing strongly interacting gauge group 
GJnew = Gs ® SU(3)c and its supcrconformal version 
GscoNF - SU(3)sc §5 SU(3)c M- 

In table (jVI| ) we have tabulated the numerical values 
of vev of H and H, the angular parameter 9, tan(0), w, 
the top mass mu and the top Yukawa coupling A^^ con- 
tributing to the parameter space of MSSM for the n — A 
level fiat directions QQQL, QuQd, QuLe and uude. It 
should be noted that appearance of large vev of Hig- 
gses as mentioned in table (jVip can easily be interpreted 
when Einstein Hilbert term appears in the total action 
of the theory at lowest order approximation [24] which 
is our present consideration. Consequently the contri- 
butions from the hard cutoff is sub-leading due to the 
soft conformal symmetry breaking. This leads to small 
top Yukawa coupling in the restricted parametric space 
of MSSM characterized by the phenomenological bound: 
43 GeV <mu < 170 GeV , 1.006 < tan{J3) < 1.025 for 
the n=4 flat directions. 
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Neglecting all the sub-leading contributions arising 
from the top Yukawa coupling in the restricted parame- 
ter space of the MSSM, the solutions of these RGE for 
n=4 level flat directions can be written as 



(37) 



gi(po) 



Here gi{^J-o), rni{p.o), A^(mo), m^ip-o) and Xpip-o) repre- 
sent the value of the gauge couplings, gaugino masses, tri- 
linear couplings, soft SUSY braking masses and Yukawa 
couplings at the characteristic scale ^o- In equation ([37]) 
we have used the following shorthand notations: 



Anii 
Ami 
Am? 



where the (3 indices 1,2,3 represent U, D, E respectively. 





i=l(U(l)Y) 


i = 2(SU(2)L) 


i = 3(SU(3)c) 


F=1(QQQL) 





3 
^ 


2 
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F=2(QuQd) 





3 
4 


8 
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F=3(QuLe) 





3 
4 


4 
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F=4(uude) 








2 
3 



TABLE III. Entries of fp matrix obtained from the solution 
of RGE 

In equation ^T^i pp and (C^)°^ are (4 x 3) and (3 x 3) 
matrices whose entries are tabulated in Table (|IIip and 
Table (jIV[) respectively. It is obvious from the RGE that 
P = 1,2 implies a = h and /3 = 3 implies a^h. 

Using the solutions of RGE along with the approxima- 
tion that the running of the gaugino masses and gauge 





i = l(U(l)v) 


i = 2(SU(2)L) 


i = 3(SU(3)c) 


/3=l(U),a=b 


26 
99 


6 


32 
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/3=2(D),a=b 


14 
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/3=3(E),a 7^ h 
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TABLE IV. Entries of (C^)"'' matrix obtained from the solu- 
tion of RGE 
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FIG. 9. Running of the ratio of the Yukawa coupling 
( A^d^o ) ) ^^'^ 1°°P I^SSM with the logarithmic 

scale logj^Q (/x). Here we have used /io = 2.6 x lO^GeV and 
/J = l([/),2(0),3(i5) Vi. 



couplings is very very small we get: 



- 



(38) 



where we have J\ — 0, J2 — 3 and J3 = 4 for i = 1, 2, 3 
and all the entries of K^'^ (3 x 3) matrix are tabulated in 
tabledV]). 

In this context the subscript '0' represents the val- 
ues of parameters at the high scale /ip. As dis- 
cussed in section III, constraining only D\ and 
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i=l(U(l)Y) 


i = 2(SU(2)L) 


i = 3(SU(3)c) 


/3=1(U) 


13 
18 


3 


8 
3 


/3=2(D) 


7 

18 


3 


8 
3 


/3=3(E) 


3 


3 






TABLE V. Entries of if'^* matrix 



is sufficient here. Egn P^ provides an extra con- 
straint relation which restricts the parameters fur- 
ther leading to more precise information in RG flow. 
For universal boundary conditions, the high scale 
is identified to be the GUT scale hgut w 3 x 
10^^ GeV, rhiifiGUT) = m2{fiGUT) = rhs^fiauT) = m, 
AEitiGur) = Au{^j,gut) = Ad{i^gut) = Aq and gi sa 
0.56, 32 ~ 0.72, gs ss 0.85. Now depending upon the 
different phenomenological situations the n = A level flat 
directions are divided into two classes. The first class 
deals with QuQd, QuLe which is lifted completely at 
n — A level. The other class which is lifted by higher 
dimensional operators deals with uude, QQQL. Most 
importantly uude, QQQL take part in the proton de- 
cay {p — 7r''e+, p TT+i^e etc.) which introduces a 
stringent constraint on the Yukawa coupling Aq at n = 4 
level. Additionally the neutrino-antineutrino oscillation 
data restricts Aq again. Then we just use RG equations 
along with these restrictions to run the coupling con- 
stants and masses to the scales as mentioned in table (jVip 
with M = 2.4 X 10^^ GeV. 



Considering all these values we obtain effectively 





« -o.osec^ 


Dl 


« -0.074C, 


Dl 


« -0.071C, 


Dl 


w -0.031C, 







-0.048 -0.168C^ 
where C — m/m^ is calculated at the GUT scale. Typ- 
ically the running based on gaugino loops alone results 
in negative values of DiVi. Positive values can be ob- 
tained when one includes the Yukawa couplings, practi- 
cally the top Yukawa, but the order of magnitude remains 
the same. The choice of fine tuned initial conditions di- 
rectly shows more fine tuning is required compared to 
other models. It is a straightforward exercise to verify 
that even if one considers all the flat directions at n = A 
level one will arrive at the potential ean. (ll3l) with same 
Co and C4. This is precisely what we have done in this 
paper. 

The results of RG flow have been demonstrated in 
figs®-®. In fig® and fig® 'dashed, 'solid' and 'dot- 
dashed' line represents U(1)y, SU(2)l and SU(3)c 
gauge group content respectively. Fig®-fig® explic- 
itly showing the behavior of the RGE flow of gaugino 
masses, soft SUSY breaking mass, trilinear couplings and 
Yukawa couplings respectively. Additionally flg®-flg® 
give consistent GUT scale unification. 



V. SUMMARY AND OUTLOOK 

In this article we have proposed a model of infla- 
tion in the framework of MSSM with new flat direc- 
tions using saddle point mechanism. We have demon- 
strated how we can construct the effective inflationary 
potential in the vicinity of the saddle point starting from 
n = A level superpotential for the flat direction content 
QQQL,QuQd,QuLe and uude for MSSM. The effec- 
tive inflaton potential around saddle point, resulting from 
the non-vanishing fourth derivative of the original po- 
tential, has then been utilized in estimating for the ob- 
servable parameters and confronting them with WMAP7 
dataset using the publicly available code GAME, which 
reveals consistency of our model with latest observations. 
We have then explored the possibility of Primordial Black 
Hole formation from the running-mass model by estimat- 
ing the mass of PBH. 

Subsequently, we have engaged ourselves in flnding out 
the effective parameter space and the constants appear- 
ing in the saddle point analysis for the MSSM inflation 
by solving the one loop RGE. It is worth mentioning that 
the RGE flow of fourth level MSSM is exactly solvable 
in this context and we hope that all the numerics can 
be tested in the LHG or any linear collider in near fu- 
ture. Gonsequently we conclude that fourth level MSSM 
inflation confronts extremely well with WMAP7 within a 
certain parameter space obtained from one loop MSSM 
RGE flow. 

A detailed survey of RG flow with two loop beta func- 
tion, inflection point inflation 26] for n = A level MSSM 
candidates, sensitivity in the neighborhood of the sad- 
dle point with the one loop corrected potential, the ef- 
fect of quantum Goleman De Luccia tunneling 27] and 
the inflationary model building of MSSM derived from 
string theory via braneworld using several compactifica- 
tion schemes remain an open issue, which may even pro- 
vide interesting signatures of MSSM infiation. We hope 
to address some of these issues in due course. 
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Flat 
direction 


Mo = (t>o 
GeV 


GeV 


ra.4,0 
GeV 


(H> 

GeV 


(H> 

GeV 


e 

in deg 


tan{e) 


w 


V 

GeV 


mu 
GeV 


GeV 


QuLe 


2.6 X 10^ 


36.967 


7.546 


0.200 X 10^^ 


0.458 X 10^*' 


45.171 


1.006 


0.994 


0.500 X 10" 


43 


1.212 X 10"" 


QuQd 


2.6 X lO'^ 


36.967 


7.546 


0.450 X 10" 


0.423 X lO^'' 


45.370 


1.013 


0.987 


0.601 X 10^^ 


170 


7.106 X 10"" 


QQQL 


1.344 X 10" 


892 X 10^ 


182 X 10^ 


0.188 X 10* 


0.124 X 10** 


45.707 


1.025 


0.975 


0.226 X 10** 


80 


4.945 X lO"** 


uude 


2.896 X 10" 


4.142 X lO'* 


845 X 10^ 


0.174 X 10^ 


0.157 X 10*^ 


45.549 


1.019 


0.981 


0.235 X 10*^ 


135 


8.047 X 10"" 



TABLE VI. MSSM parameter values obtained from RG flow for n=4 level flat directions 



[1] A. Mazumdar and J. Rocher, Phys. Rept. 497 (2011) 85. 
[2] J. D. Barrow, R. Bean and J. Magueijo, MNRAS 316 
(2000) L41. 

[3] R. AUahverdi, J. G. Bellido, K. Enqvist and A. Mazum- 
dar, Phys. Rev. Lett. 97 (2006) 191304. 

[4] A. Chatterjee and A. Mazumdar, arXiv:1103 .5758l 

[5] D. H. Lyth, JCAP 0704 (2007) 006; R. AUahverdi, K. 
Enqvist, J. G. Bellido, A. Jokinen and A. Mazumdar, 
JCAP 0706 (2007) 019. 

[6] K. Enqvist, A. Jokinen, S. Kasuya and A. Mazumdar, 
Phys. Rev. D 68 (2003) 103507. 

[7] CERN Large Hadron Collider, for more details see the 
online link: http://public.web.cem.ch/public/en/lhc 

[8] WMAP collaboration, D. N. Spergel et al., Astrophys. J. 
Suppl. 170, 377 (2007); for uptodate results on WMAP, 
see http://lambda.gsfc.nasa.gov/product/map/current 

[9] Planck collaboration, 

http://www. rssd. esa. int/index.php ?project=Planck 
some early results are also available, see, for example, P. 
A. R. Ade et.al., arXiv:1101.2022 
[10] T. Gherghetta, C. Kolda and S. P. Martin, Nucl. Phys. 

B 468 (1996) 37. 
[11] A. R. Liddle and D. H. Lyth, Cosmological Inflation 
and Large-Scale Structure, Cambridge University Press 
(2000). 

[12] CAMB, Online link: http://camb.mfo/ 
[13] A. Kosowsky and M. S. Turner, Phys. Rev. D 52 (1995) 
1739; N. Duechting, Phys. Rev. D 70 (2004) 064015. 



[14] M. Drees and E. Erfan i, JCAP 1104 (2011) 005; M. Drees 

and E. Erfani, larXi^ lllO.6052 
[15] E. J. Copeland, A. R. Liddle, J. E. Lidsey and D. Wands, 

Phys. Rev. D 58 (1998) 063508. 
[16] B. J. Carr, Astrophys. J. 201 (1975) 1. 
[17] J. Niemeyer and K. Jedamzik, Phys. Rev. Lett. 80 (1998) 

5481; J. Niemeyer and K. Jedamzik, Phys. Rev. D 59 

(1999) 124013; M. Shibata and M. Sasaki, Phys. Rev. 

D 60 (1999) 084002; J. C. Hidalgo and A. G. Polnarev, 

Phys. Rev. D 79 (2009) 044006. 
[18] H. P. Nilles, Phys. Rept. (1984) 110. 
[19] S. P. Martin, hep-ph/9709356 

[20] L. E. Ibanez and C. Lopez, CERN preprint TH-3650 
(1983). 

[21] K. Inoue, A. Kakuto and S. Takeshita, Kyushu Univ. 

preprint 83 - HE - 6 (1983). 
[22] C. Kounnas, A. B. Lahanas, D. V. Nanopoulos and M. 

Quiros, CERN preprint TH - 3651 (1983). 
[23] T. Kobayashi, H. Nakano and H. Terao, Phys. Rev. D 71 

(2005) 115009. 

[24] N. Kaloper, L. Sorbo and J. Yokoyama, Phys. Rev. D 78 
(2008) 043527. 

[25] R. Harnik, D. T. Larson, H. Murayama and M. 

Thormeier, Nucl. Phys. B 706 (2005) 372. 
[26] S. Hotchkiss, A. Mazumdar and S. Nadathur, JCAP 1106 

(2011) 002; K. Enqvist, A. Mazumdar and P. Stephens, 

larXiv: 1004.3724; R. AUahverdi, B. Duttaand Y. Santoso, 

iFXiv: 1004.2741 
[27] S. Coleman and F. De Luccia, Phys. Rev. D 21 (1980) 

3305. 



